Adsorption Trajectories and Free-Energy Separatrices for Colloidal Particles in 

Contact with a Liquid-Liquid Interface. 



Joost de Graaf/'B Marjolein Dijkstra, 1 and Rene van Roij 2, 

1 Soft Condensed Matter, Debye Institute for Nanomaterials Science, 
Utrecht University, Princetonplein 5, 3584 CC Utrecht, The Netherlands 
^Institute for Theoretical Physics, Utrecht University, 
Leuvenlaan 4, 3584 CE Utrecht, The Netherlands 
(Dated: February 18, 2010) 

We apply the recently developed triangular tessellation technique as presented in [J. de Graaf et 
al, Phys. Rev. E 80, 051405 (2009)] to calculate the free energy associated with the adsorption of 
anisotropic colloidal particles at a flat interface. From the free-energy landscape, we analyze the 
adsorption process, using a simplified version of Langevin dynamics. The present result is a first 
step to understand the time-dependent behavior of colloids near interfaces. This study shows a wide 
range of adsorption trajectories, where the emphasis lies on a strong dependence of the dynamics 
on the orientation of the colloid at initial contact with the interface. We believe that the observed 
orientational dependence in our simple model can be recovered in suitable experimental systems. 

PACS numbers: 82.70.Dd, 45.50.Dd, 68.03.Cd 



I. INTRODUCTION 

The adsorption of colloidal (nano)particles at liquid- 
liquid interfaces is not only of scientific interest, but also 
relevant for industry. The formation of two-dimensional 
structures [l|, which may be utilized in photonic 
bandgap materials, and the stabilization of Pickering 
emulsions 3], are two examples of possible applications. 
There is also an impetus to theoretically describe the 
colloidal adsorption process more accurately, in order to 
gain a deeper insight into the mechanisms at work in 
experimental systems. Especially two-dimensional (2D) 
fluid phenomena, and phase transitions in 2D fluids of 
anisotropic particles [1, [f| are of interest. The tunabil- 
ity and variety of colloidal particles currently available, 
coupled with our still limited knowledge on colloidal ad- 
sorption, leaves the study of their interfacial phenomena 
an open field. 

To better understand the complex systems which 
arise when colloidal particles are brought in contact with 
an interface, we developed a method to determine the 
free energy associated to the adsorption of a single parti- 
cle [13]. This free-energy calculation is based on similar 
surface tension arguments as proposed by Pieranski 
in his ground breaking study of colloidal adsorption phe- 
nomena. The model in our investigations encompasses 
surface and line tension, but disregards interfacial de- 
formation and electrostatic effects. Our contribution to 
this theoretical description is the triangular tessellation 
scheme [ltj, which allows us to efficiently determine the 
adsorption free energy of an anisotropic colloid. In this 
paper we apply our technique to perform an initial in- 
vestigation into the dynamics that occur when a colloid 
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attaches at the interface and relaxes to its preferred con- 
figuration. 

The dynamics we study in our model is a simplified ver- 
sion of Langevin dynamics [ill fl2j . From the adsorption 
free energy, which acts as a potential energy on the col- 
loid, a vector field of the adsorption force is determined. 
This vector field is then studied by flow lines. For conve- 
nience, the anisotropy of the particle is only considered 
in the free energy and not in the friction tensor, which 
for interfacial systems would also be position dependent. 
Our choice allows us to examine general trends and gross 
features, and to showcase our numerical method, whilst 
at the same time give results which resemble those one 
would expect for actual physical systems. We consider 
the time scale associated with the behavior predicted by 
our model, and conclude that there are parameters for 
which the dynamics are sufficiently slow to allow these 
effects to be observed in experimental systems. 

In this paper we focus on three systems to point out the 
complexities that anisotropy can induce in the adsorption 
process. As a reference system we discuss the adsorption 
of an ellipsoidal particle for several conveniently chosen 
system parameters. Subsequently, we proceed to describe 
the same system for a cylindrical particle to show the ef- 
fects of shape. The free-energy landscape for cylinders 
is more complicated, leading to metastable adsorption 
configurations with a large stability domain in the ad- 
sorption region. Finally, we study a cylindrical colloid 
with aspect ratio 1. The occurrence of interfacial ad- 
sorption here is strongly dependent on the way in which 
the particle makes initial contact with the interface. For 
some initial configurations there will be no adsorption at 
the interface, depending on the value of the three-phase 
contact angle of the particle and the two liquids. 

In conclusion, we apply our triangular tessellation 
scheme to determine the adsorption behavior of non- 
convex particles via flow-line dynamics. We show that 
for the systems considered the adsorption mechanism is 
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far more complex than was previously believed. The ef- 
fects of including interfacial deformation, and anisotropic 
friction, certainly merit further investigation, but go be- 
yond the scope of our initial investigation. 

II. METHOD 
A. Theoretical Considerations 

We consider the adsorption of a solid uniaxial convex 
colloid at a planar oil-water interface separating two ho- 
mogeneous half spaces of oil and water. The interfacial 
normal coincides with the z axis of our coordinate frame, 
see Fig. [1] The depth z of the interface is measured with 
respect to the center of the particle. We locate the origin 
of our system at the particle's center. Medium 1 (Mi) is 
the half space above and medium 2 (Ms) is the half space 
below the interface. The polar angle <p € [0,7r/2] gives 
the angle between the colloid's rotational symmetry axis 
and interfacial normal. For simplicity, capillary effects 
due to the presence of the colloid at the interface are 
neglected. A simple calculation [13] shows that gravity 
does not play a significant role in interfacial adsorption 
for most colloidal systems, and it is hence ignored. 

The adsorption free energy of the colloid contains 
terms depending on three surface areas with correspond- 
ing surface tensions: (i) the surface area of the colloid 
above the interface Si, (ii) the colloid's surface area be- 
low the interface S2, and (iii) the surface area excluded 
from the interface by the presence of the colloid S12 ■ The 
contact line, of length L, where the three phases meet, 
also contributes to the adsorption free energy. The de- 
pendence of these quantities on z and <f> is implicit for 
notational convenience. We denote the total surface area 
of the colloid by S, therefore S\,Sz € [0, S], under the 
constraint S = Si + S2 ■ 

Following Pieranksi [l[ , the adsorption free energy can 
be written as 

F{z,<f>) = 712 [(Si - S) cos 9 - S12] + tL, (1) 

where 712 is the M1-M2 surface tension, and r is the line 
tension. This free energy, Eq. (JTJ, is defined with respect 
to a reference point, namely F = when the colloid is 
completely immersed in M\. The contact angle 6 is in- 
troduced via Young's equation 712 cos6> = 7i c — 72c fl3j - 
with 7i c and 72c the Mi-colloid and the M2-colloid sur- 
face tension respectively. F(z* , <p) is made dimensionless 
and scale invariant by dividing the free energy by 712 S 
(712 7^ 0) and writing z — z*\Ja? + 2b 2 . Here a is the 
rotational symmetry semiaxis, b the perpendicular semi- 
axis, and m = a/b the aspect ratio. We thus obtain 

f(z*,<f>) = = cos9 ^ - !) - r i2 + r% (2) 

712-3 

where r\ = S\/S and ri2 = S12 / S are surface area ratios, 

r* ^ (3) 
712V0 
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FIG. 1: Two representations of an ellipsoidal colloid adsorbed 
at a flat interface, located at depth z measured from the cen- 
ter of the colloid (0,0). The xz view (a) shows the two me- 
dia, their dividing interface, and the polar angle <f> which the 
colloid's rotational symmetry axis makes with the interfacial 
normal. The surface area of the colloid above the interface is 
denoted by Si, with 7i c the Mi-colloid surface tension, and 
the surface area of the colloid below the interface is denoted 
by S2, with 72c the M2-colloid surface tension. The colloid 
excludes an area 5*i2 from the interface, the region enclosed 
by dashed curve in the xy view (b), of which the surface ten- 
sion is given by 712. The length L of this dashed curve is 
the contact line length and the corresponding line tension is 
denoted by r. The solid curve in the xy view indicates the 
colloid's outline. 



is the dimensionless line tension, and / = Lj\[S is a di- 
mensionless contact line length. An extended adsorption 
free-energy expression, which can handle nonconvex pat- 
terned colloids, is given in Ref. [Toj [ . 

The adsorption configuration is the location of the free- 
energy minimum in Eq. ([2]) and is denoted by (z* d , </> a d)- 
The adsorption free energy for this configuration is given 
by /ad = /(^ad'^ad)- For multiple (metastable) min- 
ima, the various (z* d ,</> a d) are labeled with a subscript 
i = 1,2,..., where the deepest minimum is given the 
lowest index. We denote the positive value of z*, for 
which the interface just touches the top of the particle 
by z^ et (4>) for given 0; as it is also the detachment posi- 
tion. Note that this value is always positive, whenever we 
consider — zi et {<p) the interface just touches the bottom 
of the particle, because of the symmetry of the prob- 
lem. When z* < -Za et (<f>) or z* > z^ et ((f>) the colloid 
is completely immersed in M\ and M2 respectively. For 
a given <f>, f(z*,(j)) has a minimum as a function of 2*, 
the location of which is denoted by zj^ in ((/)) and the cor- 
responding free energy by /min^)- It is possible that a 
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single equi-0-curve f(z*,<p) has multiple minima. This 
in turn leads to multiple z* nin (</>), which are labeled with 
indices as above. Analogously, <j) m i-a{z*) gives the loca- 
tion of the minimum in f(z*, </>) as a function of <fi for a 
fixed z* , the value of which we denote by / m in(z*)- Here 
multiple minimum curves can exist as well. Often we ab- 
breviate .2d ct 0), 0min(^*), and z* lin 0) to z*[ et , <f> min and 
z min: respectively. 



B. Numerical Approximation Scheme 

Deriving analytic expressions for the dependence of Si , 
S2, S\2, and Lonz' and 4> is highly nontrivial in gen- 
eral, if not impossible [13]. To analyze colloids adsorbed 
at an interface we have developed the numerical tech- 
nique called triangular tessellation, which we briefly out- 
line here, see Fig. [5] Full details are given in Ref. [To| . 
including the application to more complex colloid shapes 
than discussed in this paper. 

Using a parametrization of the colloid at the interface, 
the surface and shape of the particle are approximated 
by a polyhedron, the faces of which are triangles, see 
Fig. (Hb). This approximation can be improved upon 
by reducing the overall triangle size, thus achieving bet- 
ter correspondence between the actual particle shape and 
the approximating polyhedron. The fineness of the tes- 
sellation, which is often noncquidistant, is indicated by 
annxm vertex notation [Io| . The value nm is the num- 
ber of vertices in the mesh, each of which is a member 
of 6 adjacent triangles. Modeling a 2D or 3D object in 
this way is a well-known computer science method and it 
has been successfull y appli ed to various surface tension 
problems in physics jl4l - [l6j . The area of each of the tri- 
angles composing the polyhedron can be determined by 
a simple cross product. Summing over these areas results 
in an approximated value for the total surface area, we 
denote this value by S. 

When the tessellated object is intersected by a plane, 
some triangles of the approximating polyhedron are in- 
tersected by it, see Fig. HJb). These triangles are subdi- 
vided, Fig. (He), into at most three subtriangles, which 
are not intersected [13] • After this subdivision, we deter- 
mine which triangles lie above, given by the set Af, and 
below the plane, given by the set A4., see Fig.[2ji. Using 
these sets it is possible to determine the approximated 
surface areas Si = X^^T,«> ano - ^2 = zCi^^i, where 
A*,i is the surface area of the z-th triangle in the set A*. 
It follows that Si + S 2 = S. 

From the polyhedral approximation the points are ex- 
tracted where the plane intersects ribs of the triangles 
that tessellate the colloid, see Fig.[5Jc,f). These are used 
to obtain the approximate surface area excluded from the 
interface by the colloid, say Si 2, and the approximate 
contact line length, say L. The points of intersection 
form a polygon in the plane, Fig. the area of which 
is determined using a trapezium rule type integration or 




FIG. 2: A sketch of the triangular tessellation scheme in a xz 
plane projection (a-d) and in a xy plane projection (e,f). A 
colloid at a planar interface (a) is approximated by a polyhe- 
dron of which the faces are triangles (b). The surface of the 
colloid is thus effectively tessellated by triangles, which can 
be used to determine the approximate total surface area S. 
Triangles which are intersected by the interface are divided 
into at most three subtriangles (c). Upon this subdivision, 
both the area above and below the interface is approximated 
more accurately (d). Using the sets of triangles in graph (d), 
it is possible to determine Si and 52, which approximate the 
value of 5*i and S2 respectively. The boundary of the area cut 
out from the interface, the dashed line in (e) , is approximated 
by a polygon (f ) . This polygon is formed by the points where 
the interface intersects triangles of the tessellation (b-c). The 
surface area S12 and contact line length L can be determined 
from (f). 



a modified version of Green's integral theorem [10( . With 
either method, it is important to take steps to minimize 
numerical uncertainty [10|. 



C. Adsorption Trajectories and Separatrices 

Applying the triangular tessellation technique to solve 
Eq. J2) as a function of z* and </> yields a free-energy 
landscape (on a grid). This landscape is studied in vari- 
ous ways, using z*- or 0-sections [13], and by ^in(^) an d 
curves with respective corresponding /,„;„. It is 
also possible to examine the negative gradient at each 
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point of the adsorption free-energy landscape 
F{z\<t>*) = -Vf(z 4 ,f); 

with z* and <f>* unit vectors, 7r</>* = <ft, and f*{z*,<f>*) = 
/(z*, 0); thus obtaining a vector field of adsorption force 
T . Note that with this choice for </>* and /* the domain 
is bounded by (z*,<f>*) G [-1,1] x [0,0.5]. It should be 
emphasized that the adsorption free energy acts as a po- 
tential energy for the colloid, which can be differentiated 
with respect to its macroscopic coordinates, obtaining a 
force; the microscopic coordinates of the fluid have been 
integrated out to yield the tension terms. 

The vector field is studied by examining four features: 
(i) flow lines, (ii) minima and maxima, (hi) saddle points, 
and (iv) separatrices. We will show that for certain par- 
ticle species there are multiple minima in the adsorption 
free-energy landscape, also see Ref. [lCf. Each of these 
minima is surrounded by a region to which that par- 
ticular minimum is attractive. That is to say, all flow 
lines which originate from points in this region reach 
that minimum when t — > oo, with t time. A flow line, 
r](t) = (z*(t), <j)*(t)), is defined to follow the path of steep- 
est descent from its starting point at rj(t — 0). It is a 
solution to the differential equation 

m = ^-=nv(t)), (5) 

where the dot denotes the time derivative and the so- 
lution is fixed by imposing the initial point rj(t = 0). 
We refer to these flow lines as adsorption trajectories to 
stress their relation with the physical path followed by 
a particle when adsorbing to the interface. The dividing 
line between two attractive regions is referred to as the 
separatrix. This separatrix can contain (local) maxima 
and saddle points. 

D. Timescale of Colloid Motion 

It is desirable to review the concept of time as intro- 
duced by solving the differential equation f]{t) — F{rj{t)). 
The variable t in this equation is mathematically a 
parametrization parameter. This parameter is however 
related to physical time, since it is relative to the size of 
the force exerted at a point 77(f), but not scaled according 
to a friction tensor. In the Langevin equation [ll|: E3 f° r 
a sphere in a homogeneous medium, the time to travel 
over a short distance Al is equal to M£AZ/|A|, where M 
is the mass, the Brownian time, and A is the force. 
We have a similar situation for our parametrization, save 
the mass and friction coefficient. The 1/|.F| is the re- 
duced time that it takes to move along a flow line through 
a specific point. Note that it diverges for extrema in the 
free-energy landscape, which is in accordance with the 
particle being stationary there. We analyze the relation 



between the reduced time f and the actual time t r in 
more detail in the "Discussion" section. 

The flow line for a specific initial point is determined 
by linear differential solver which employs steepest de- 
scent. In this scheme the time parameter is determined 
as follows. The time in which step i is taken, Afj, is 
defined as the ratio of the step size AZj and the force at 
the initial point of the step, i.e., Af, = AZ,/|V/j|. We 
thus obtain a series of Af,-. The time required to arrive 
at point n along the flow line, say f(«), from the initial 
point is then given by 

n-1 

t(n) = 5^A*i. (6) 
i=i 

In the limit of infinitesimally small step size, the solution 
to the vector field with the initial condition rj(t = 0) and 
the proper time dependence are obtained. 

III. RESULTS 

All data have been produced using the triangular tes- 
sellation technique 10] on 500 2 -vertex grids heteroge- 
neously mapped to the object under consideration. To 
determine the landscape, the value of f*(z* ,</>*) is cal- 
culated on 250 by 250 nonequidistant grid points in 
z* G [-Zd 0t , Zd 0t ] and 0* e [0,0.5]. The data on this 
mesh is then interpolated with a 3rd order interpolation 
scheme to yield the full landscape. The accuracy of the 
triangular tessellation has been verified by a semianalytic 
method [lOj. In this case, the technique has been proven 
to have a relative uncertainty of less than 10~ 4 per data 
point. The flow lines are determined by means of a linear 
differential solver, and contain over 1000 steps per line. 
The step size is reduced until subsequent reductions re- 
sult in a relative difference between the curves of less 
than 10~ 3 . This uncertainty is determined by summing 
over all points, the difference in position between two suc- 
cessive reductions and dividing by the total path length. 
The differential solver is considered to have converged on 
the solution of the vector field when the difference is less 
than 10 -3 . For some extreme cases we are forced to work 
with less precision, as we will indicate in the text. 

A. Adsorption of Ellipsoidal Particles 

Figure[3Ja) shows adsorption trajectories 77(f) for an el- 
lipsoidal particle with aspect ratio m = 6 for cos 9 = —0.5 
and r* = -0.1. In Fig. ^b) the ±z^(tj>), the < in (0), 
the </> m in(z*), and the attractor A(z*,(f>) curves for this 
landscape are given. The location of the minimum is in- 
dicated with a dot. The term attractor is introduced here 
to describe a feature in the vector field of adsorption force 
to which the flow lines are attracted, as can be seen in 
Fig. EJa). Mathematical analysis shows that for points 
on the attractor one of the eigenvectors of the Hessian 
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FIG. 3: (Color online) Graph (a) shows adsorption trajecto- 
ries r/(t) for an ellipsoidal particle with aspect ratio m = 6 
for cos# = —0.5 and r* = —0.1. The arrow heads indicate 
the direction of colloid motion through the free-energy land- 
scape, the dot gives the location of the minimum, i.e., the ad- 
sorption configuration, and the thick gray lines the ±z*j e t(</0 
curves. The symbols M\ and M2 indicate that the colloid is 
completely immersed in the respective media, whenever its z* 
value is to the right or left of the ±z dot (0) curve. Graph (b) 
shows the ±z det {4>) curves in gray (thick, solid), the z^ irL {(j)) 
curve in red (thick, dots), the <?imm(z*) curve in green (thick, 
dash-dot), and the attractor A(z*,(j>) in blue (thin, solid). 
The location of the minimum is again given by a dot. Note 
that the £mm(</>) an< A the <j)min{z*) curves in (b) indeed in- 
tersect with points on the adsorption trajectories where the 
tangent is vertical and horizontal respectively, when (a) and 
(b) are superimposed. 



matrix [(VV T )/] is in the direction of the gradient, with 
a positive sign, and the other is perpendicular to it, with 
a negative sign. However, the gradient is nonzero for 
points on the attractor, which are therefore not saddle 
points. Note that the attractor lies between the z^ in (0) 
and </> m in(-z*) curves, for which one of the components of 
the gradient vector is zero. 



By analyzing all points on the landscape, we conclude 
that the entire region between the ±z*j ct (</>) curves with 
4> G [0,7r/2] is attracted to the single minimum with 
(^ad'^ad) ~ (—0.0368, 0.57r). This was to be expected 
on the basis of our observations in Refs. [IOJ. The flow 
lines in Fig.^a) give a rather abstract picture of the col- 
loidal motion through the interface. The main point of 
introducing this ellipsoidal system is to have a basis for 
comparison, when we study cylinders in the next section. 

To illustrate the behavior of the particle along one of 
the adsorption trajectories in Fig.^a), we have included 
Fig. 21 which shows the time dependent movement. Fig- 
ure IDJa) shows several snapshots of the motion along the 
flow line with r](t = 0) w (0.930, 0.0965tt) in FigH Fig- 
ure U{b) shows the location of the snapshots on the flow 
line. Here it should be stressed that the minimum is 
reached only at t 00. The final snapshot shows a con- 
figuration, which is reached at a finite time, very close to 
the minimum. 

From Fig. Ufa) it becomes clear that the movement 
of the colloid through the interface, when it is close to 
vertical, is quite slow when compared to the rotational 
part of the movement before it reaches its equilibrium 
configuration. The final part of the movement, reaching 
the equilibrium position, takes infinitely long however. 
Slowing down at the end occurs only when the colloid is 
very close to its adsorption configuration. The difference 
in speed between the rotational and vertical motion parts 
of the adsorption is caused by the difference in free-energy 
decay between points to 7 and 8 to 15 in Fig @|b) 
respectively. In part it is caused by our implicit choice of 
the ratio between friction coefficients for the translational 
and rotational part of the movement. We return to this 
in the "Discussion" section. 

The results presented here are for one contact angle 
cos 9 and one reduced line tension r* only. Further anal- 
ysis and our observations in Ref. [Io| show that qualita- 
tively similar results are obtained for contact angles with 
-1 < cos# < and r* < 0.0. For r* > the formation 
of adsorption barriers, as discussed in Ref. [10j, hinder 
adsorption to the interface form either medium and con- 
sequently change the adsorption free-energy landscape 
sufficiently to alter the adsorption trajectories qualita- 
tively. 

Additional representations of the movement along the 
adsorption trajectories shown in Fig. |3ta) are given in 
the Appendix. The motion studied in Fig. [4] will also be 
revisited. 



B. Trajectories for Cylindrical Colloids 

Figure [5] shows the adsorption trajectories and asso- 
ciated curves for a cylindrical colloid with m = 6 when 
cos0 = —0.5 and t* = —0.1. Despite the close similarity 
of these parameters to those in section A for ellipsoids, 
the corresponding free-energy landscape for cylinders has 
two minima instead of one: at (z* d , </> a d) ~ (—0.079, 0.57r) 
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FIG. 4: Graph (a) shows several snapshots of the motion of an 
ellipsoidal colloid with m — 6, cos# = —0.5, and r* = —0.1 
through the interface along a flow line from Fig [3] An xz 
view of the colloid is represented by a black silhouette and 
the interface is indicated by a thin black line. The time, at 
which a snapshot is taken for each frame, is given in the top 
left corner. The numbers at the bottom right of the interfacial 
line, correspond to the numbers in graph (b). The first frame 
also shows the location of the two media Mi and M2 . Graph 
(b) shows the ±2;| et (0) curves, thick gray line; the adsorption 
free-energy minimum, black dot (overlapped by point 15) ; the 
flow line for which the snapshots are taken, thin black curve; 
and the (z* , <j)) for which the snapshots in graph (a) are taken, 
numbered black squares. Again the location of the media is 
indicated by the symbols Mi and M2. 



with / ad « -0.397 and at (z* d ,<?(> ad ) w (-0.962, O.Ott) 
with / a d ~ —0.090. For these parameters the mini- 
mum with </> a( j = 0.5-7T is the absolute minimum, not un- 
like the ellipsoidal case of section A, the minimum with 
<^ad = 0.07T is a metastable minimum. It can be shown 
theoretically that the latter constitutes an configuration 
where one of the caps of the cylinder is flush with the 
interface. Thereby a surface is excluded from the inter- 
face, lowering the free energy, whilst the rest of the colloid 
is immersed in the preferred medium (Mi in this case). 
There is only contact with M2 on the excluded surface. 
This configuration can be metastable, because excluding 
surface area from the interface strongly lowers the free 
energy. 

The kinks in the adsorption trajectories of Fig. Eta) 
are directly related to ridges in the free-energy landscape 
caused by the sharp corners of the cylinder. These ridges 
are indicated by the z^ sc (4>) curves in Fig. [S^b), where 
"dsc" refers to the discontinuity that occurs in the gra- 
dient vector field. Note that similar to the z^ ct curves 
these are also symmetric in z* =0, hence we will use the 
± notation to distinguis between the different branches. 
In the bottom left corner of Fig. EKb) the three curves 
nearly coincide, but it can be shown that from left to 
right they are m in(z*), z^ n {(j>), and -z^ sc ((f>). From the 
figure it is not clearly apparent that the curve 4> m in{z*) is 
discontinuous at (z* , </>) (0.315, O.Ott). There is a small 
gap between the 'vertical' and 'horizontal' branches. The 
kink in the vertical branchand the starting point of the 
horizontal segment lie above each other. For zj^ in (0) the 
ridges in the potential landscape also induce a kink. The 
presence of these ridges can strongly influence the behav- 
ior of the flow lines for the cylindrical colloid's adsorption, 
as we will see. 

In Fig. Efc), we show the separatrix S(z*,<f>) between 
the two minima and the three attractors Ai{z*,(f), i = 
1,2,3, which are present in this free-energy landscape. 
There is a saddle point on the separatrix, where it meets 
with two of the attractors, A\ and Ai- The separatrix 
forms the division between the regions to which the re- 
spective minima are attractive. Remarkably the size of 
the (z*,(f>) domain to which the metastable minimum 
is attractive, is still substantial for an aspect ratio as 
high as m = 6. For a particle adsorbing to the interface 
from Mi, on the left-hand side of Fig. 03c), any colloid 
which touches the interface with <f> < 0.2247T will adsorb 
to the metastable minimum. For attachment from M2 
we find that the metastable minimum is attractive for 
4> < 0.0147T. Therefore, the colloid has only a very small 
window to adsorb to the metastable configuration from 
the energetically unfavorable medium. However, adsorp- 
tion to the metastable configuration from the preferred 
medium is quite likely, because almost half of the orien- 
tations will lead to this configuration. 

The appearance of a tertiary attractor, see Fig. [5jc), 
is somewhat surprising. The first and second attractor 
are merely the split form of a feature similar to the at- 
tractor of the ellipsoid from section A. This "main at- 
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FIG. 5: (Color online) Properties of the free-energy landscape 
for a cylinder with aspect ratio m — 6 for cos 9 = —0.5 and 
r* = —0.1. Graph (a) shows several adsorption trajectories 
rj(t) (thin black curves); the location of free-energy minima 
(dots), of the saddle point (diamond), and of media Mi and 
A/2- Thick gray curves indicate ±Zdct(0)- Graph (b) displays 
the vector field discontinuities, #dsc(0)> m purple (thick, dots). 
The Zmin(^) curve is represented in red (thick, dashed) and 
0min(2*) in green (thick, dash-dot). Graph (c) shows the 
attractors, Ai(z* ,</)), and the separatrix, S(z*,(f>). 



tractor" is split, because of the two minima, causing the 
separatrix to intersect it. The tertiary attractor, leading 
to the absolute minimum of the free-energy landscape, 
is caused by a subtle interplay between the right-most 
ridge, z*[ sc (0), and the 4> m in(z*) curve. Note that the at- 
tractor terminates at exactly the same point as where the 
0min(z*) curve has a kink. Also note that the presence 
of the tertiary attractor strongly influences the behavior 
of the adsorption trajectories around it, see Fig. [5fa). 

For this cylindrical colloid with m = 6, we again expect 
qualitatively similar results for r* < 0; for r* > adsorp- 
tion barriers are found. However, since the depth of the 
minima is strongly dependent on the value of the contact 
angle cos 9, also see Ref. [Tc| for this depth when m = 4, 
the position of the separatrix will vary significantly with 
cos 9. We come back to this in the next section, when we 
discuss cylindrical particles with aspect ratio m = 1. 

The adsorption trajectories in Fig. [5ja) that end up 
in the stable minimum, are similar to those of an el- 
lipsoid, except around the secondary attractor, see the 
Appendix for examples. In Fig. [5] we consider col- 
loid movement along one flow line which leads to the 
metastable minimum. Figure H][a) reproduces several 
snapshots of the motion along the flow line with rj(t — 
0) « (0.979, 0.014tt), see also Fig [5] Figure [^b) shows 
the location of the snapshots on the flow line. The final 
snapshot gives a configuration very close to the mini- 
mum, which is reached only at infinite time. It should 
be noted, that the motion between frames 12 and 15 (see 
Fig. IS]) could be determined only with a relative uncer- 
tainty of R3 10~ 2 . This is due to the close proximity of the 
</>min(-2 : *), z^ in (<j>), an d — z dsc('/ > ) curves, as can be seen in 
Fig.[5j Consequently, the time dependence of this part of 
the motion has significant uncertainty, which we estimate 
to be at most 15%. 

Similar to the ellipsoid in Fig. @}a) the cylinder first 
moves almost vertically through the interface (frame 1 — 
8), before it tilts slightly (frame 9 — 11). However, when 
it touches the interface with one point of the edge of one 
of the end caps, see frame 12, it does not continue to 
tilt in the same direction. At this point it is energetically 
favorable to move into the metastable minimum by tilting 
back to the vertical orientation (frame 13 — 15). Finally 
the colloid comes to rest with one end cap flush with the 
interface at t — > oo. Tilting back when the edge of the 
end cap makes contact with the interface, frame 12, is an 
indication that the discontinuity ridges in the landscape 
act as a dynamical "barrier" . 

We have thus shown that if a particle can exclude a 
relatively large area from the interface by one of its ends, 
the dynamics and possible adsorption configurations are 
strongly influenced. Both an ellipsoid and a spherocylin- 
der do not have a secondary minimum and therefore can- 
not exhibit the vertical adsorption of single particles. 
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FIG. 6: Graph (a) shows several snapshots of a cylindrical 
colloid with m = 6, cos 6 — —0.5, and r* = —0.1 adsorb- 
ing to the metastable minimum. An xz view of the colloid 
is represented by a black silhouette and the interface is indi- 
cated by a thin black line. The time is given in the top left 
corner (first two rows) and the bottom left corner (last two 
rows). The numbers at the bottom right of the interfacial 
line, correspond to the numbers in graph (b). The first frame 
indicates the location of the two media. Graph (b) shows the 
±Zdet(4>) curves, thick gray line; the adsorption free-energy 
minimum, black dots (one is overlapped by point 15); the 
saddle point, black diamond; the flow line for which the snap- 
shots are taken, thin black curve; and the (z*,<p) for which 
the snapshots in graph (a) are taken, numbered black squares. 



C. Separatrices and Special Configurations 

In this section we study the free-energy landscape of a 
short cylindrical colloid with aspect ratio m = 1, r* = 
and several cos0. Figure [TJ a) shows the minima, sad- 
dle points, and separatrices for five values of the contact 
angle, cos(0) = 0.0, -0.1, -0.2, -0.4, and -0.6. In 
Fig. \T[b) these three properties of the free-energy land- 
scape are given for cos(0) = —0.7, —0.72, —0.74, and 
—0.78. The inset shows the separatrices for a small re- 
gion of Fig. |7fb). Finally, Fig. [7[c) indicates the loca- 
tion of the minima, saddle points, and separatrices for 
cos(0) = -0.8, -0.9, and -0.95. The inset is included 
to prove that the various minima in the top left corner 
of the graph are all separate points, as opposed to the 
minima in the bottom left corner of the graph, which all 
have the same z* and 4> value. 

For cos = there is a single minimum at (z* d , 4>&d) ~ 
(0.000, 0.27l7r) and for cos0 = —1 there is no minimum. 
For all other values of cos one of the minima is located at 
(z* d , <^ad) ~ (—0.571, 0.0-71"), which is exactly when one of 
the cylinder caps is flush with the interface. The location 
of the other minimum changes with the value of cos 0. 
One of these minima is stable and the other is metastable. 
The minimum with (z* d ,^ a( j) ~ (— 0.571, O.Ott) is stable 
when cos0 < —0.23. The labelling order of the minima 
with a( j = 0.57T in the tops of Figs.JTJb) and (c) indicates 
the way in which these appear in the graph. As already 
mentioned, the inset of Fig.[7|c) shows that these minima 
are indeed distinct. 

With decreasing cos0 (cos0 tends towards —1) the 
separatrices, denoted by the black lines, move from the 
left lower corner towards the center of the region en- 
closed by the ±z det (</>) curves, the adsorption region. 
The location of the secondary minimum shifts closer to 
(z* d ,0 ac j) ~ (—0.571, 0.571"). For a certain value of cos0 
the separatrix splits into two pieces. This is shown in 
Fig-EJb), where the separatrices for cos0 = —0.7, —0.72, 
—0.74, and —0.78 are given (labels 5, 6, 7, and 8, respec- 
tively). Further study indicates that for cos0 ps —0.705 
there is a transition between the single and double sep- 
aratrix regime. From Figs. [TJa) and (b) it is apparent 
that for decreasing cos0 the position of the saddle point, 
denoted by a diamond, moves closer to the —z^ et (4>) 
curve. Around cos = —0.705 the saddle point lies on the 
— z dot (0) curve and is degenerate. For a normal saddle 
point there are two attractive and two repulsive direc- 
tions, whereas here there are two repulsive and only one 
attractive direction. 

For lower values of cos0 the saddle point transforms 
into two degenerate saddle points on the boundary of the 
adsorption region. Both of these have only one attractive 
and one repulsive direction. The inset in Fig. EJb) clearly 
shows that the separatrices are disjoined over the entire 
adsorption region. This was to be expected, since the sep- 
aratrices are flow lines for a vector field. FigureEJc) gives 
a clearer picture of the distance between the two sepa- 
ratrices for cos0 < —0.705. To determine these curves a 
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FIG. 7: Properties of the free-energy landscape of a cylin- 
der with aspect ratio m = 1 for r* = and several contact 
angles, labeled as follows in (a) cos# = 0.0 (0), —0.1 (1), 
-0.2 (2), -0.4 (3), -0.6 (4); in (b) cos (9 = -0.7 (5), -0.72 (6), 
-0.74 (7), -0.78 (8); in (c) cos6> = -0.8 (9), -0.9 (10), 
—0.95 (11). Minima are given by thick dots, saddle points 
by diamonds, separatrices by black lines, and the ±«| et (</>) 
curves by thick gray curves. The inset in (b) shows an en- 
largement of a piece of graph (b). The location of some of the 
minima in graphs (b) and (c) are shown in the inset in (c). 



tolerance of 10~ 2 was used for the flow line convergence 
algorithm. When cos 8 J, — 1 the area between the two 
separatrices tends to encompass the entire adsorption re- 
gion. The most interesting feature of this inter-separatrix 
domain is that a flow line starting in it will flow towards 
a point on the — z det ((/>) curve, between the two saddle 
points. Along such a flow line the value of f(z*,<p) is 
always positive, yet monotonically decreasing. 

The effect of an inter-separatrix domain for certain val- 
ues of the contact angle is that the adsorption of colloids 
can be strongly dependent on the initial configuration. 
The existence of adsorption minima does not necessar- 
ily imply that adsorption will take place, even if thermal 
fluctuations are ignored. However, the situation studied 
here, a cylindrical colloid with aspect ratio m = 1 and 
extreme values of the contact angle, is not representative 
for most colloidal systems. To illustrate the differences 
between the various "adsorption" possibilities for a cylin- 
der with m = 1 and cosO = —0.95, we have included 
Fig. HI which shows snapshots of the movement of the 
colloid for three different adsorption trajectories though 
the free-energy landscape. Figure [SJd) gives the location 
of the various snapshots in Figs|5fa)-(c). For clarity the 
rotational symmetry axis is given by a white line on the 
black silhouette of the colloid. 

The first trajectory, Fig. [5][a), shows the colloid ad- 
sorbing to the metastable minimum at (.z* d ,</>ad) ~ 
(—0.567, 0.57r). In its final configuration, the colloid only 
barely penetrates the interface and its rotationally sym- 
metry axis is perpendicular to the interfacial normal. 
Figure |5Jc) shows the adsorption to the primary mini- 
mum at (z* d ,0 a d) = (—0.571, 0.071"). Here the final con- 
figuration is when one of the cylinder caps is flush with 
the interface, as discussed earlier. Note that the rota- 
tional symmetry axis is indeed parallel to the interfacial 
normal. The intermediate series of snapshots, Fig. Efb), 
shows the colloid along a flow line in the inter-separatrix 
region. The colloid moves through the interface and slows 
down when it approaches the — z*l ei (4>) curve, since the 
gradient tends to zero here. The final configuration in- 
cluded shows the cylinder with only a single attachment 
point on the interface, it is in essence detached. The three 
trajectories shown in Fig.[8]have also been included in the 
movie files referred to in Appendix. 



IV. DISCUSSION 

In the previous we have seen that the dynamics of col- 
loid motion through the interface can be very rich. In 
this section we analyze the relation between the pro- 
posed flow-line dynamics and the dynamics associated 
to experimental systems. We also discuss the elements 
present in real systems, which we have not incorporated 
in our simple model and consider how they might affect 
the obtained results. 
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FIG. 8: (Color online) Graphs (a-c) show snapshots of the motion of a cylindrical colloid with m = 1, cos 9 = —0.95, and 
r* = 0.0 through the interface along a flow line of the corresponding free-energy landscape. An xz view of the colloid is 
represented by a black silhouette and the interface is indicated by a thin black line. The white line on the silhouette is used 
to indicate the rotational symmetry axis of the colloid. The time is given in the top left corner (first row) and the bottom left 
corner (second row). The numbers at the bottom right of the interfacial line, correspond to the numbers in graph (d), the graph 
letter is given near the "7" on each flow line. For clarity graph (d) shows the three respective flow lines r/(t) corresponding to 
the snapshots in graphs (a-c). The flow lines are indicated by a thin black curve, the separatrices S(z*,(f>) [also see Fig. [TJc)] 
by a thin yellow curve, the ±z^ et (<j>) by a thick gray curve, and the minima by a thick dot. In all graphs the location of the 
two media is given by the symbols Mi and M2. 



A. Estimate of the Adsorption Time 

For a dispersed particle undergoing a force, the dynam- 
ics are governed by a solution to the complete Langevin 
Equation [ll|: E3 ■ When we neglect inertia, the random 
force term, and limit ourselves to studying the z and 4> 
components, the equations of motion (EOMs) become 

9^) = _d p 
at r oz 

= -^F(z(t r ),cf>(t r )), (8) 

where t r is the "real" unreduced time, F is our unsealed 
adsorption free energy, and a is the long semiaxis of the 



particle. The prefactors A and /1 represent the trans- 
lational and rotational friction coefficients respectively. 
These coefficients should account for the size a and the 
shape of the particle, as well as, the positional z and ori- 
entational <f> dependencies of the system. This includes 
the difference in viscosity between the two media, which 
influences the friction force, the way in which friction 
changes in the z direction according to the colloid's </>, 
and vice versa. We have not included these dependencies 
in the notation used here, because of the simplifications 
we apply to Eqs. © and §5$ in the following. 

Let us assume that A and \i are independent of z and <fi; 
therefore, only changing with scaling and particle type; 
we return to the validity of this later. It can then be 
easily shown that the EOMs of Eqs. © and © allow 
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us to regain the flow-line dynamics of Eq. ([5]), 77(f) 
J"(r7(f)). We require that F{z,(j>) = ^ X2 Sf*{z*, (/>*), z 

Va 2 + 2b 2 Z*, <f> = TT(f>*, t r 



nt, where k is a time scale 



based on the system parameters, and 



fi(a) 



2b 2 



a 2 + 2b 2 

7125 



A(a); 
A(a), 



(9) 
(10) 



where the friction coefficients are now coupled by impos- 
ing the flow-line criterion. 

For a sphere we set A(a) = 6nriQa, with 770 the viscosity 
of the densest medium. In this way we regain the ex- 
pected translational friction [ItJ when the particle is free 
to move in that medium. Note that for this choice, we 
do not recover the expected rotational friction coefficient 
87T77oa 3 [l8|, [ll| in front of the time derivative d<j)/dt r . 
Instead we find 18?7oa 3 /7r, which differs by a factor of 
9/47T 2 ss 0.23. In retrospect, it would have been more 
appropriate to rescale the system such that the flow-line 
dynamics give the expected results for a sphere. How- 
ever, we investigated the adsorption trajectories with this 
constraint and found qualitatively similar results. 

Before discussing the accuracy of the simplifications 
used in our research, let us gauge the time t r required for 
a colloid to attach to the interface and relax to its final 
position. We assume A(a) = 67r?7oa, with 770 the largest of 
the two viscosities. For the movements studied in Figs.|H 
|6l and [8] the duration of the adsorption process (the part 
in which significant changes occur) is of the order f rs 5. 
According to Eq. (fT0|). the physical time is estimated at 



6717700- 



26 2 



7125 



-t. 



(11) 



Upon using this approximation we arrive at the follow- 
ing result, see Table HI Here we have considered an ellipse 
with aspect ratios m — 1 and m = 6, sizes a = 125 nm, 
a = 500 nm, and a = 2.5 /im, which are typical col- 
loidal length scales. For the viscosity we use 770 = 1.5 
Pas (glycerol), 770 = 1.56 • 10 -2 Pas (cyclohexylchlo- 
ride), and 1.0 -10~ 3 Pas (H2O) [2(j, which are commonly 
used solvents. Surface tension values between two typ- 
ical liquids are of the order 712 = 10~ 2 Nm -1 [21| or 
lower than 712 = 10~ 5 Nm -1 with the addition of sur- 
factants [H[23||. 

Note that it follows from Table U that length of the 
adsorption process can vary significantly with the choice 
of system parameters. For values of the interfacial ten- 
sion which are in the order 10 2 Nm -1 the adsorption 
process typically takes tens to hundreds of microseconds. 
However, upon lowering the interfacial tension, e.g., by 
adding surfactants, the process slows down significantly, 
in some cases taking upwards of a second. This is an 
exciting prospect, since for many systems the time scale 
is thus accessible to experimental techniques. It can also 
be shown that the particular shape of the particle, e.g, 
ellipsoid, cylinder, or spherocylinder, does not influence 



these results significantly, since (a 2 + 2b 2 ) /S is virtually 
the same for these particle types. 

The limiting factor in an experiment to determine the 
adsorption behavior is the maximum operating frequency 
of the camera and not the optical elements used for imag- 
ing. For modern cameras operating frequencies of 10.000 
Hz or more are obtainable [24]]. However, we should bare 
in mind that our time t r is based on a rather crude esti- 
mate and that there can be a substantial deviation from 
the prediction from this value for an actual experimen- 
tal system. Nevertheless, a discrepancy of a decade or 
two will still place many of the above mentioned systems 
within the experimentally observable range. 



Limitations of the Model and Experimental 
Considerations 



There are several key points of criticism which can be 
identified when trying to compare the dynamics predicted 
by our flow lines in the reduced system to real- world dy- 
namics. These points should be taken into account in 
follow-up studies, but go beyond the scope of an initial 
investigation. 

Neglecting the random thermal force component is an 
acceptable simplification, if we only want to study the av- 
erage movement, i.e., the trajectory averaged over many 
adsorption events. Some care should be taken though, 
since it is clear from, for instance, Fig. [TJb) that ad- 
sorption trajectories which have almost identical starting 
points can diverge from each other quite rapidly. This 
may prove problematic in establishing an average for an 
experimental system, especially when the magnitude of 
the free-energy landscape features becomes in the order 
of a k#T. Neglecting the inertia is also acceptable, since 
most colloidal systems are strongly overdamped. Even 
under high potential differences — a 1 /im colloidal par- 
ticle can typically experience an interfacial adsorption 
potential in the order of 10 4 -keT [l| — colloids cannot 
be easily forced into the ballistic regime of motion. 

Another point of concern is neglecting anisotropic ef- 
fects in friction force. As we have seen in the previous, the 
friction coefficient in part imposes the time scale on the 
system. For anisotropic particles, which have a friction 
coefficient tensor, this time scale changes according to the 
direction of motion and the position of the particle. How- 
ever, the difference in friction coefficient between motion 
along the long semiaxis of an ellipsoid and perpendicular 
to this semiaxis has been shown to be no greater than 
a factor of two [12, [25l - [27j . For the rotational friction, 
the change is likely to be more substantial than a factor 
of two for strongly oblong particles. Suppose we assume 
proportionality with the rotational radius to the third 
power, as is the case for a sphere Then, depending 
on the axis of rotation, the friction coefficient varies be- 
tween /i oc a 3 and /1 cx 6 3 . Even if the proportionality is 
less extreme than in our simple estimate, the differences 
can easily be very significant. 
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TABLE I: The physical time t r required to complete the motion of the colloid through the interface on the basis of our model. 
We assume t ~ 5 and study two aspect ratios m — 1 and m — 6 for ellipsoids. Several values of the viscosity rjo are considered. 
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The difference in friction between the two media bor- 
dering the interface, can be several orders of magnitude, 
see for instance Ref. [2(| for the viscosity values of com- 
monly used solvents. Such a difference in friction coeffi- 
cient will play a role in describing a particular system ac- 
curately. We assumed that the timescale of the dynamics 
are dominated by the medium with the highest friction 
coefficient, or, equivalently, viscosity. This assumption 
holds whenever there is substantial penetration of the col- 
loid in that medium. A more accurate approach requires 
us to weight the friction coefficient by the surface areas 
of the particle which are in the respective media. Strong 
differences in viscosity of the respective media most likely 
induce a resistance for the colloid to move into the more 
viscous medium. This can substantially alter the quali- 
tative behavior of the flow lines. 

For experimental systems the solvent is sometimes 
density-matched to the particle to eliminate the effects 
of gravity. Matching is however more complicated to 
achieve for a three component system, such as the parti- 
cle near a liquid-liquid interface. To account for gravity 
in our model, whenever, we require the buoyancy mass 
of the particle. This mass can be determined by vol- 
ume integration over the parts which penetrate the re- 
spective media, assuming a homogeneous density distri- 
bution. Such volumetric integrations can be performed, 
using a 3D analogy of our 2D surface integration scheme 
based on triangular tessellation. That is to say, some 
kind of spatial tessellation using polyhedra. Neverthe- 
less, we have shown in Ref. [10| that for colloidal systems 
the effects of gravity will be negligible, because of the low 
Bond number for such particle sizes. 

In many physical system, the particle and the inter- 
face are often charged by self dissociation of the surface 
molecules on the colloid and of those which compose the 
media [28]. For a like-charged interface there can be a 
substantial charge repulsion barrier which needs to be 
overcome before adsorption can take place, if this bar- 
rier can be crossed at all. Including electrostatic effects 
into the model will substantially change the appearance 
of the free-energy landscape. Finally, we should stress 
that including interfacial deformation [2^, [3(| into our 
model should qualitatively change the results, whenever 
the capillary effects are dominant. Together with electro- 



static effects, the results of deformation cannot be pre- 
dicted at this time. 

In conclusion, we have identified the time which we 
expect the adsorption process to take for several experi- 
mentally feasible systems, based on our model. The re- 
sults are encouraging as this time regime seems to be ac- 
cessible to current experimental techniques, even though 
there are still quite a few caveats to be considered. 

V. CONCLUDING REMARKS 

We have employed the triangular tessellation technique 
introduced in Ref. (lpj to determine the adsorption free- 
energy landscape of anisotropic colloids at a flat inter- 
face. This free energy is composed of surface and line 
tension contributions, but does not take into account in- 
terfacial deformation and electrostatics. We analyzed the 
obtained free-energy landscapes by means of a vector field 
of adsorption force and its associated flow lines. These 
flow lines are calculated using a linear steepest descent 
method and they are parametrized in units of reduced 
time. The link between these flow line dynamics and 
Langevin dynamics and the validity of our approxima- 
tion are examined. For typical viscosities and surface 
tensions we predict that timescales to complete the ad- 
sorption process are in the microsecond to second regime. 
Such timescales are accessible for observation in modern 
experimental set-ups. 

Within our model, we re-establish that there is a strong 
dependence of the adsorption free-energy landscape on 
the shape of the colloidal particle. This in turn leads to 
a wide range of adsorption phenomena. For ellipsoidal 
particles we find a single equilibrium adsorption config- 
uration with relatively simple dynamics. However, for 
cylindrical particles there can be two adsorption config- 
urations, depending on the value of the contact angle. 
The presence of two minima in the free-energy landscape 
leads to a separatrix, which forms the divide between 
the regions to which the respective minima are attrac- 
tive. The metastable minimum can have a large domain 
in the adsorption region to which it is attractive. This 
is significant, since it suggests that adsorption in unex- 
pected configurations should be easier than previously 
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believed. The implications of this large domain size on 
the experimentally observed buckling and flipping tran- 
sitions [3, Q merits further investigation. 

We have also shown that for short cylinders the free- 
energy landscape allows for another type of colloid in- 
teraction with the interface, other than adsorption. De- 
pending on the value of the contact angle, there can be 
a domain in the adsorption region, for which colloids do 
not attach to the interface, but simply move through it 
unhindered. This is of particular interest, since it shows, 
within the confines of our model, that the presence of 
stable and metastable adsorption configurations is not a 
sufficient criterion to guarantee particle adsorption. To 
the best of our knowledge, this is the first time that this 
phenomenon has been observed theoretically. 

Our simple model to calculate the free energy does not 
include interfacial deformation due to capillary, electro- 
static, or gravitational forces; nor does it take into ac- 
count inertia, anisotropy in the friction tensor, and the 
Brownian random force to determine the colloid dynam- 
ics. We are however optimistic that the phenomenol- 
ogy described here can be recovered in a more elaborate 
model. The method presented here thus forms the first 
stepping stone towards obtaining better correspondence 
between experiments and their theoretical description, 
and thus brings us to a closer understanding of colloidal 
adsorption phenomena. 
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Appendix A: Colloid Motion Though the Interface 

In this Appendix we describe the content of the EPAPS 
files [3l| for this paper. The EPAPS files give movies for 



the time dependent movement along a few adsorption 
trajectories representative of the colloid's motion through 
the free energy landscapes of Figs. 02 [5] and [51 



The data repository contains three folders labeled as 
follows: the particle type appears first, followed by the 
aspect ratio m, followed by the value of the contact an- 
gle cos0, and finally the line tension value r* is given. 
The lower case "m" appearing in front of the numbers, 
is to indicate that it is a negative value, where the "m" 
stands for minus. The individual movie files are labeled 
similarly: an abbreviation of the type, followed by the 
value of the aspect ratio, and finally the coordinates of 
r](t = 0) in the order z*(t = 0), <f>*(t = 0). The file where 
the value of the aspect ratio is preceded by "A" give the 
motion of the colloid along the attractor. 



Each movie file consists of a panel which shows two 
frames side by side. The left frame gives the ±z*[ et lines 
in gray as well as the flow line in black. The final con- 
figuration is indicated by a black dot and the configu- 
ration at time t by a red dot. The time appears in the 
top left corner of the right frame. A thin line divides 
this frame in two and serves to indicate the location 
of the interface. A silhouette similar to those used in 
Figs. SJa), IU a ) and[8ja-c) gives the xz projection of the 
colloid. Again a white line is added to the silhouette of 
the cylindrical colloid with m — 1 to show the colloids 
orientation. 



A few of the movie files show slight jittering of colloid 
along its adsorption trajectory. Especially flow lines to 
the metastable minimum of a cylinder with m = 6, also 
see the main text. These adsorption trajectories suffered 
from slightly higher levels of numerical uncertainty than 
the others, since many features in the landscape are close 
together in this case. 
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